Pre-equilibrium dilepton production from an anisotropic quark-gluon plasma 
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We calculate leading-order dilepton yields from a quark-gluon plasma which has a time-dependent 
anisotropy in momentum space. Such anisotropics can arise during the earliest stages of quark-gluon 
plasma evolution due to the rapid longitudinal expansion of the created matter. Two phenomeno- 
logical models for the proper time dependence of the parton hard momentum scale, Phard, and the 
plasma anisotropy parameter, ^, are constructed which describe the transition of the plasma from 
its initial non-equilibrium state to an isotropic thermalized state. The first model constructed inter- 
polates between 1-1-1 dimensional free streaming at early times and 1-1-1 dimensional ideal hydrody- 
namical expansion at late times. In the second model we include the effect of coUisional broadening 
of the parton distribution functions in the early-time pre-equilibrium stage of plasma evolution. 
We find for both cases that for fixed initial conditions high-energy dilepton production is enhanced 
by pre-equilibrium emission. When the models are constrained to fixed final pion multiplicity the 
dependence of the resulting spectra on the assumed plasma isotropization time is reduced. Using 
our most realistic coUisionally-broadened model we find that high-transverse momentum dilepton 
production would be enhanced by at most 40% at the Relativistic Heavy Ion Collider and 50% 
at CERN Large Hadron Collider if one assumes an isotropization/thermalization time of 2 fm/c. 
Given sufficiently precise experimental data this enhancement could be used to determine the plasma 
isotropization time experimentally. 

PACS numbers: 11.15Bt, 04.25.Nx, ll.lOWx, 12.38Mh 



I. INTRODUCTION 

With the ongoing ultrarelativistic heavy ion collision experiments at RHIC and future experiments planned at LHC 
the goal is to produce a deconfined plasma of quarks and gluons (QGP). This new state of matter is expected to 
be formed once the temperature of nuclear matter exceeds the critical temperature, Tc ~ 200 MeV. However, many 
properties of the expected QGP are still poorly understood. One of the most difficult problems is the determination of 
the isotropization and thermalization time of the QGP, Tiso and rtherm, respectively.^ At RHIC energies it was found 
that, for ^2 GeV, the plasma elliptic flow, W2, was well described by ideal hydrodynamical models. Based on early 
studies [J |3j 13 ideal hydrodynamical fits to elliptic flow data indicated that the matter can be modeled as a nearly- 
perfect fluid starting at extremely early times after the collision, riso ~ 0.6 fm/c However, recent hydrodynamical 
studies Q which include the effect of all 2nd-order transport coefficients consistent with conformal symmetry have 
shown that these initial estimates for the isotropization/thermalization time of the plasma have a sizable uncertainty 
due to poor knowledge of the proper initial conditions (CGC versus Glauber), details of plasma hadronization and 
subsequent hadronic cascade, etc. As a result, it now seems that isotropization times up to Tjso ^ 2 fm/c are not 
ruled out by RHIC data and in order to further constrain this time additional theoretical and experimental input will 
be required. 

As mentioned above there are uncertainties introduced in hydrodynamic modeling of the quark-gluon plasma 
(QGP) due primarily to the dependence of the results on the assumed initial conditions: initial energy density, spatial 
profile, fiow velocity, etc. In order to remove this uncertainty one would like to find observables which are sensitive 
to the earliest times after the collision and are relatively unaffected by the later stages of plasma evolution. One 
obvious candidate to consider is high-energy dilepton production since dileptons couple only electromagnetically to 
the plasma and therefore, in their high-energy spectra, carry information about plasma initial conditions. To this end 
here we calculate the dependence of leading-order high-energy dilepton production on the assumed isotropization and 
thermalization time of the QGP using two simple models for early-time QGP evolution. 

To begin the discussion we introduce two proper time scales: (1) the parton formation time, tq, which is the time 



For simplicity from here forward we will assume that these two time scales are the same, r^hcim = 'Hso) so that the system achieves 
isotropization and thermalization at the same proper time given by TigQ. 
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after which one can treat the partons generated from the nuclear colhsion by a distribution of on-shell partons; and 
(2) the isotropization time, Tiso, which is the time when the system becomes isotropic in momentum-space. At RHIC 
energies tq ~ 0.3 fm/c and at LHC energies tq ^0.1 fm/c. Immediately after the collision, the partons are produced 
from the incoming colliding nuclei at r = tq, at which time the partonic momentum distributions can be assumed to 
be isotropic.^ The subsequent rapid longitudinal expansion of the matter (along the beam line) causes it to become 
much colder in the longitudinal direction than in the transverse direction Longitudinal cooling occurs because 
initially the longitudinal expansion rate is larger than the parton interaction rate and, as a result, a local momentum- 
space anisotropy is induced with (pf) ^ (p^) in the local rest frame. If the system is to return to an isotropic 
state it is necessary that at some later time the interaction rate overcomes the expansion rate with the system finally 
isotropizing and remaining isotropic for r > tisq. Once isotropy is achieved (and maintained by parton interactions) 
the use of hydrodynamic simulations can be justified. We are therefore critically interested in the properties of the 
plasma around the isotropization time as these provide the relevant initial conditions for subsequent hydrodynamic 
evolution. 

The study of anisotropic plasmas has received much interest recently due to the fact that a quark-gluoiiplasma 
which has a local momentum-space anisotropy, 2 (iOr. ) ^ (2t)jJs subject to the chromo-Weibel instability 0, H 0, [13, 
M, [3 El, Q d, M, [13, M, M, H HH, m, M, MM, III, [23. The chromo-Welbel instability causes rapid growth 
of soft gauge fields which preferentially work to restore the isotropy of the quark-gluon plasma on time scales much 
shorter than the coUisional time scale. However, most of the theoretical and numerical developments in describing 
the time-evolution of a QGP subject to the chromo-Weibel instability have been restricted to asymptotic energies 
at which perturbative resummations can be applied and hence the presence of the instability-driven isotropization 
at RHIC and LHC energies is not yet proven. In addition, numerical studies of the chromo-Weibel instability in an 
1-dimensionally expanding system show that there is a time delay before the effects of plasma instabilities become 
important to the system's dynamics 0, [l^]. Future work will address these issues but until they become available 
there is a substantial amount of theoretical uncertainty in the QGP isotropization time, tq < Tiso ^ 3 fm/c. 

In the absence of a precise physical framework for describing the thermalization of the quark-gluon plasma and 
the associated time scales, one possible way to proceed is by studying the dependence of observables sensitive to 
the earliest times after the collision on the assumed plasma isotropization time by constructing simple space-time 
models. As mentioned above, one candidate observable is electromagnetic radiation such as high-energy photon 
and dilepton production'^ since these particles interact only electromagnetically and can escape from the strongly 
interacting medium created after the collision unhindered. Hence, they are perfect probes for studying the early-time 
dynamics of the system. In the case of high-energy medium photon production it is difficult for experimentalists to 
subtract the large backgrounds coming from tt" decays from other sources of photons, making it hard to measure a 
clean high-energy medium photon production signal.^ In the case of high-energy dileptons, the experimental situation 
is dramatically improved and it then becomes a question of making the necessary theoretical predictions to see how 
large the effect of a possible anisotropic pre-equilibrium phase would be on dilepton production. 

Phenomenological studies of the production of high-energy dileptons have shown that there are several important 
dilepton sources and it's necessary to include each of these depending on the kinematic region. For dilepton pair 
transverse momentum or invariant mass greater than 1 GeV the most important sources of dilepton pairs are: charm 
quark decays, initial state Drell-Yan scatterings, jet-conversion, jet-fragmentation, and mediimr (thermal) production. 
For an extensive discussion of the various sources of the dilepton production in a heavy-ion collision, we refer the 
reader to [H, and references therein. Note that in most previous phenomenological studies of photon and dilepton 
production it has been assumed that the system "instantaneously" thermalizes with Tiso = tq and hence is locally 
isotropic throughout its evolution. 

In this paper, we extend our previous work 30] and concentrate on the impact of momentum-space anisotropics on 
the leading order medium dilepton production at large invariant mass and transverse momentum. We propose two 
simple phenomenological models for the time dependence of the plasma momentum-space anisotropy, £, = {Pr) I (Pl) ~ 
1, and hard momentum scale, phard- In the first model we interpolate between 1-1-1 dimensional longitudinal free- 
streaming and 1-1-1 dimensional ideal hydrodynamic expansion. In the second model we incorporate the effect of 
momentum-space broadening due to hard-hard elastic scatterings in the pre-equilibrium dynamics.^ In both models 
we introduce two parameters, Tiso and 7, with 7 setting the width of the transition from early-time pre-equilibrium 



^ See Ref. Q for an alternative view of the early times after the initial collision wherein the authors find that the distribution may be 
prolate for longer than assumed here. We postpone the study of the possibility of early-time prolate distributions to future work. 

^ Hereafter, high-energy dileptons will refer to lepton pairs with pair transverse momentum (pt) ox invariant mass (M) greater than 1 
GeV. 

* For a discussion of the effect of possible momentum-anisotropies on high-energy photon production see Ref. fsS] . 

^ We will assume that this elastic scattering rate is regulated in the infrared by an isotropic screening mass for simplicity. 
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FIG. 1: Sketch of the time dependence of (left) the plasma anisotropy parameter, ^ — ^{pt)/{p'l) ^ 1) B.nd (right) the 
normalized energy density, £ /£o, in the three different cases considered here. Note that in the right panel we only sketch the 
case of fixed initial conditions. We also adjust our models to constrain them to fixed final particle multiplicity. 



dynamics to late-time equilibrated dynamics. In the limit Tiso — > tq our models reduce to ideal 1+1 dimensional 
hydrodynamical expansion and in the opposite limit Tiso — > oo correspond to two different types of non-equilibrium 
plasma evolution: free streaming expansion or coUisionally-broadened expansion. 

In Fig. [1] we sketch the proper time dependence of the plasma anisotropy parameter, ^, and energy density in 
the three cases considered below: free-streaming followed by ideal hydrodynamic expansion, coUisionally broadened 
expansion followed by ideal hydrodynamic expansion, and "instantaneous" thermalization (ideal hydro throughout 
evolution). As can be seen from this figure for both the free streaming and coUisionally-broadened cases with fixed 
initial conditions the energy density is always greater than that obtained by a system undergoing ideal hydrodynamic 
evolution and therefore one expects that, for fixed initial conditions, non-equilibrium effects can significantly enhance 
dilepton production. In order to calculate just how much the signal is enhanced requires detailed calculations which 
we present below. As our first result we will show that for fixed initial conditions the addition of a pre-equilibrium 
phase at times t < Tiso enhances high-energy dilepton production by up to an order of magnitude. 

We then study the effect of constraining our space-time models so that the initial conditions vary in order to 
guarantee that the final pion multiplicity is independent of the assumed isotropization/thermalization time. We show 
that allowing the initial conditions to vary in this way reduces the final effect of the possible anisotropic pre-equilibrium 
phase. In order to quantify the effect we introduce a ratio called the "dilepton enhancement" , (j), which is the ratio 
of the dilepton yield obtained when one assumes a finite isotropization/thermalization time to that obtained when 
one assumes that the plasma "instantaneously" thermalizes at the formation time, Tiso = tq. We show that, in our 
most realistic coUisionally-broadened model, the dilepton enhancement can be up to 40% at RHIC energies and 50% 
at LHC energies if one assumes a isotropization/thermalization time of 2 fm/c. We also show that as one varies 
the assumed isotropization/thermalization time that the dilepton enhancement, 0, has a non-trivial dependence on 
pair transverse momentum which could allow experimental determination of Tjso if the experimental medium dilepton 
yields are obtained with high enough precision. 

The work is organized as follows: In Sec. |TT] we calculate the dilepton production rate at leading order using an 
anisotropic phase space distribution. In Sec. IIIII we construct models which interpolate between anisotropic and 
isotropic plasmas. In Sec. IIVI we present our final results on the dependence of dilepton production on the assumed 
plasma isotropization time and compare with other relevant sources of high-energy dileptons. Finally, we present our 
conclusions and give an outlook in the Sec. |Vl 



II. DILEPTON RATE FROM KINETIC THEORY 



From relativistic kinetic theory, the dilepton production rate diV'^' /d^Xd^P = dB}^^ /d^P (i.e. the number of 
dileptons produced per space-time volume and four dimensional momentum-space volume) at leading order in the 
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electromagnetic coupling, a, is given by [31U32I. I33l|: 



Upi)fMv,,a[f 6^^\P~Pi -V2), (2.1) 



d^P J (2^)3 (27r)3 ^"^-^^■'"^-^^ "•"'"99 

where is the phase space distribution function of the medium quarks (anti-quarks), Vqq is the relative velocity 

/+; 



between quark and anti-quark and cr' is the total cross section 



where m; is the lepton mass and M is the center-of-mass energy. Since we will be considering high-energy dilepton 
pairs with center-of-mass energies much greater than the dilepton mass we can safely ignore the finite dilepton mass 
corrections and use simply cr^g' = Aira'^ /3M'^ . In addition, to very good approximation we can assume that the 
distribution function of quarks and anti-quarks is the same, fq = fq- 

In this work we will assume azimuthal symmetry of the matter in momentum-space so that the anisotropic 
quark/anti-quark phase distributions can be obtained from an arbitrary isotropic phase space distribution by squeezing 
(^ > 0) or stretching (^ < 0) along one direction in the momentum space, i.e. 

/,,,-(p,e,Phard) = /^f|(v/pMl(p^,Phard) , (2.3) 

where phard is the hard momentum scale, n is the direction of the anisotropy and ^ > — 1 is a parameter that reflects 
the strength and type of anisotropy. In general, Phard is related to the average momentum in the partonic distribution 
function. In isotropic equilibrium, where ^=0, Phard can be identified with the plasma temperature T. To give another 
specific example, in the case of 1-1-1 dimensional free-streaming discussed in Sec. IIII A 21 phard is given by the initial 
"temperature" Tq. 

For general ^ we split the delta function in Eq. (|2.ip such that we can perform the p2 integration: 

dR'^'' 



d^P 



5q2 




727r5 J 
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f d^Pi 


72tt^ J 


K Fj 
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P2 



/9(Pl,Phard,0 /g(P2,Phard,<^) (j'"*' (P - Pi - P2) 



/<? (Pl , Phard , /? (P - Pi > Phard , C) '^(-B - E^^ - Ep,^ ) 

Choosing spherical coordinates with the anisotropy vector n defining the z axis, we can write: 



(2.4) 



P2=P-P1 



Pi = Pi (sin COS , sin sin , COS 0p J , 

P = P(sin 6'p cos 0p, sin 6'p sin (/)p, cos 6'p) . (2.5) 
It is then possible to reexpress the remaining delta function as: 

2 

(5(i?-i?p, -i?pJ = 2(i?-pi)x-^/'e(x)5]<5((/..-0pJ , (2.6) 

i 

with X = 4p2p2 Qp sin^ — (2pi(£^ — P cos 6'p cos J — Af^)^. The angles (j)i are defined as the solutions to 
the following transcendental equation: 

2pi {E-Pcos9pCosep,)-AP 

cos 0i - = — — — ^ . 2.7 

2Ppx smC/p smC/pi 

We point out that there are two solutions to Eq. (|2.7p when x > 0. After these substitutions and expanding out the 
phase space integrals, we obtain: 



dR ' 5a 



d4p 18 



J dicOsOpJ J ^Pl fq ^y^pf (1 +CC0S2 0pJ,Phard 



xfqlJCE- pi)2 + e(pi COS^p, - PcOS0p)2,Phard ) , (2.8) 
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FIG. 2: The differential dilepton rate as a function of transverse momentum (left) and invariant mass (right). For the invariant 
mass dependence (left) we fixed pt= 3 GeV and for the transverse momentum dependence (right) we fixed M=3 GeV. In both 
cases Phard = l GeV and rapidity y—Q. 



with 

_ M2 

~ 2{E ~ Pcos{ep±9p,)) ■ *- ' 

Note that when ^ = 0, the limit of isotropic dilepton production is recovered trivially. Also note that as ^ increases we 
expect the differential dilepton rate to decrease since for fixed phaid the increasing oblateness of the parton distribution 
functions causes the effective parton density to decrease: 



ard J 



/ •^'^ (V^P^+ ^(p-n)2, Phard) , 

n{i = 0,Phard) _ Phard ^3 10) 



In order to evaluate the anisotropic dilepton rate it is necessary to perform the remaining two integrations in 
Eq. (|2.8p numerically. In Fig. [2] we plot the resulting differential dilepton rate as a function of transverse momentum 
and invariant mass for ^ e {0,1,10,100}. One can see the effect of increasing ^ for fixed Phardj namely that the 
dilepton production rate decreases due, primarily, to the density effect mentioned above. 

Knowing the rate, however, is not enough to make a phenomenological prediction for the expected dilepton yields. 
For this one must include the space-time dependence of phard and ^ and then integrate over the space-time volume 

^y=^Rlj d^PT I y_ ^rdrd, , (2.11a) 

where Rt — 1.2 A^^^ fm is the radius of the nucleus in the transverse plane. These expressions are evaluated in 
the center-of-mass (CM) frame while the differential dilepton rate is calculated for the local rest frame (LR) of the 
emitting region. Then, the dilepton pair energy has to be understood as E^r = px cosh {y — rj) in the differential 
dilepton rate dR^nn/d'^P . Additionally, in Eqs. (j2.1ip we have assumed that there is only longitudinal expansion of 
the system. Since at early times the transverse expansion is small compared to the longitudinal expansion, one can 
ignore it. Some studies have suggested that the influence of the transverse expansion on the space-time evolution 
becomes phenomenologically important around 2.7 fm/c [s^l, therefore, our approximation is valid for describing the 
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early-time behaviour we are interested in. Substituting Eq. (|2.8p into Eqs. (|2.1ip we obtain the dilepton spectrum 
including the effect of a time-dependent momentum anisotropy. 

Note that we have not included the next-to-leading order (NLO) corrections to the dilepton rate due to the complex- 
ity of these contributions for finite ^. These affect dilepton production for isotropic systems for E /T < 1 [ssl. [sgI [371 [sst . 
In the regions of phase space where there are large NLO corrections, we will apply if-factors to our results as indi- 
cated. These i^T-factors are determined by taking the ratio between NLO and LO calculation for an isotropic plasma, 
therefore, in this work we are implicitly assuming that the X-factor will be the same for an anisotropic plasma. 

III. SPACE-TIME MODELS 

In this section we present two new models for dimensional non-equilibrium time-evolution of the QGP and 
review the cases of 1-1-1 dimensional free-streaming and 1-1-1 dimensional hydrodynamic expansion. In all cases 
considered below the number density will obey ?t,(t) oc (tq/t) in its asymptotic regions.^ This results in all cases from 
the assumption that the total particle number is fixed while the size of the box containing the plasma is expanding 
at the speed of light in the longitudinal direction (Id expansion). 

In each case below we will be required to specify a proper time dependence of the hard-momentum scale, phard, and 
anisotropy parameter, ^, which is consistent with this scaling for r <C Tiso and r ^ Tiso- Before proceeding, however, 
it is useful to note some general relations. Firstly we remind the reader that the plasma anisotropy parameter is 
related to the average longitudinal and transverse momentum of the plasma partons via the relation 

Therefore, we can immediately see that for an isotropic plasma that ^ = 0, and for an oblate plasma which has 
> 2{pI) that C > 0. 

Secondly we note that given any anisotropic phase space distribution of the form specified in Eq. (j2.3p the local 
energy density can be factorized via a change of variables to give 

r £;3p 

f (Phard,0 = J P /iso ( + C (P • ) , phard ) , (3.2) 

= fo(Phard)'7^(f) , 

where £q is the initial local energy density deposited in the medium at tq and 

We note that hmj_»o7?.(^) — 1 and hm^^oo'^(0 — 

A. Asymptotic Limits of the Anisotropic Phase Space Distribution 

Before presenting our proposed interpolating models we review previous calculations for the free streaming and 
hydrodynamic expansion cases [3l] . [39l . [40j and show how to determine our anisotropic phase space distribution 
function parameters, phard and ^, in these two cases. 

1. 1+1 Dimensional Ideal Hydrodynamical Expansion Limit 

We first consider the limiting case that Tiso = tq so that the plasma is assumed to be "instantaneously" thermal 
and isotropic and undergoes ideal 1+1 dimensional hydrodynamical expansion throughout its evolution. In ideal 



The interpolating models will only obey this relation outside of a region of order 7^^7-150 around Tiso where the transition between 
different types of expansion takes place. In the transition region n will increase due to non-equilibrium effects as wc discuss later in the 
text. 
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hydrodynamical evolution using the boost-invariant 1+1 Bjorken model (4l| we can identify phard with the temperature 
and the anisotropy parameter vanishes by assumption, ^ = 0. Due to the fact that ^ = the distribution function 
for highly relativistic particles will depend only on the ratio between the energy and temperature, /hydro (p, 2^) — 
f{E/T{T)) with E = {p^+vlf^. In this case the number density, hard scale (temperature), energy density, and 
anisotropy parameter obey the following 



n{T) ^ no (^^^ , 

Phard(T) = T{t) = To 



Sir) 







(3.4a) 

(3.4b) 

(3.4c) 
(3.4d) 



In order to obtain an analytic result for the differential dilepton rate which is applicable at high energies one can 
approximate the quark and anti-quark Fermi-Dirac distributions by Boltzmann distributions and integrate Eq. (j2.8p 
analytically. In this case it is also possible to perform the necessary integration of the rate over the plasma space-time 
evolution analytically [sil, [3§| to obtain: 



dN, 



hydro 



dydM"^ 
dNt'" 



5a^ 1 
6^ M4 



hydro 



dRPd^prdy 




(3.5a) 



(3.5b) 



■ + z^)Ko{z) +AzKi{z) (4-|-z2)^ G{z) = z^ {8 + z^) K^iz) and mr = y^AP+p^. As a check of 



where H{z) = . 

our numerics we have verified that numerical integration of our dilepton rate given in Eq. 
Eqs. p. lip reproduces this analytic result in the limit Tjso — tq and ^ = 0. 



over space-time via 



2. 1+1 Dimensional Free Streaming Limit 

As another limiting case we can assume instead that our 1+1 dimensional expanding plasma is non-interacting. If 
this were true then the system would simply undergo 1+1 dimensional free-streaming expansion (3ll. [40j. Since, in 
this case, the system would never become truly thermal or isotropic this corresponds to taking the opposite limit from 
the one we took in the previous subsection, namely we will now take the limit Tiso 00. 

In the free streaming case, the distribution function is a solution of the collisionless Boltzmann equation 

P ■ dxfi.s.{p,x) = , (3.6) 

where the subscript f.s. indicates that this is the free-steaming solution. In this work we will also assume that the 
distribution function is isotropic at the formation time, r = tq. 



/f.s.(p,a;) 



^AMIA], (3.7) 



Phard 



where px is the transverse momentum, p^ is the longitudinal momentum and phard is the hard momentum scale at 
Tq. The typical hard momentum scale of particles undergoing 1+1 dimensional free streaming expansion is constant 
in time. In the case of indefinite free-streaming expansion the system never reaches thermal equilibrium and so the 
system strictly cannot have a temperature associated with it; however, since our assumed distribution function is 
isotropic at r = tq, we can identify the initial "temperature" of the system, Tq; with the hard momentum scale phard 
when comparing hydrodynamic and free streaming expansion. 

Eq. (|3.6p has a family of solutions which are boost invariant along the z (beam) axis 

h...{p,x) ^ f{pT,PLt ^ Ez) . (3.8) 

Therefore, the functional dependence of the distribution function for the free streaming case is of the form 

= /(vSIS^EllZM) . ,3...) 
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This distribution function can be simphfied if we change to co-moving coordinates: 

PL = Pt sinhy , E = px coshy , (3.10a) 

z — T sinh?7 , t = t cosh 77 , (3.10b) 

where, as usual, y is the momentum-space rapidity, r is the proper time, and 77 is the space-time rapidity. In terms 
of these variables one obtains 



/f.s.(p,x)=/(^^^l + ;^sinh2(y-,y)J . (3.11) 

Note that in the case of indefinite free-streaming at late times the quark and anti-quark longitudinal momentum are 
highly red-shifted reducing late time emission of high-energy dilepton pairs. 

As written in Eq. (j3.ip the anisotropy parameter is related with the average transverse and longitudinal momenta 
of the partons. The average momentum-squared values appearing there are defined in the standard way: 

2 ^ _ /rf^pp|,^/(p,x) 

Using the 1-t-l dimensional free streaming distribution given in Eq. p. lip and transforming to co-moving coordinates 
defined in (|3.10p so that d^p p\ cosh y dpr dy we obtain 



(pf)f.s. cx2^o^ (3.13a) 

„2 \ rp2 ^0 



(pi)f.s. (xT^^. (3.13b) 



Inserting these expressions into the general expression for ^ given in Eq. p.ip one obtains ^/.^.(r) = t^/tq — 1. 
With this in hand we can also determine proper time dependence of the energy density in the free-streaming case 
by substituting this expression for ^ into Eq. (|3.2p . £i.s.{t) = ^oT^i^f.s. {'''))■ early times one must use the full 
expression given by Eq. (j3.2p : however, at late times one can expand this result to obtain ff.s.(T) oc tq/t as expected 
for a 1-1-1 free streaming plasma 42]. 

Summarizing, one finds in the 1+1 free streaming case that in the limit t ^ tq: 

n(r) = no (^) , (3.14a) 

Phard(T) = ?3hard(T = 0) = Tq , (3.14b) 

£{t) = £o (^) , (3.14c) 



r . 

e(T) = 4-1- (3-14d) 

With the distribution function given by Eq. p. lip , the dilepton spectrum can be calculated. As a function of the 
invariant mass M , one obtains 



dydA'P ^ 727r3^^^^^° / '^^^ ^{t/to) 

X [{xiX2f - (xia;2COsh(yi - yz) - l/2)2]-i/2 



In the last expression, the integration is over all Xi from to +00 and over all yi from —00 to 00 subject to the 
constraint 

1 „ 1 

^ 2xiX2 ^ 



cosh(yi - y2) + 1 cosh(yi - ya) 
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As a function of the transverse momentum, px, the dilepton production using free streaming case we obtain^ 
""^'-[YoKK-mJ +^ (— smh2;2-xsmhyg j j 

with 

M 

a;± = —7 -— ^ r — ; r- . (3-17) 

2 (ttit cosh (2/1 - 7/2) ± Pt) 

We have verified that using the expressions listed in Eq. (|3.14p our direct numerical integration of the rate given in 
Eq. (|2.8p over space-time via Eqs. (|2.1ip reproduces this analytic result in the free-streaming limit. 

We note in closing that as a solution of the coUisionless (non-interacting) Boltzmann equation, the free-streaming 
case can be taken as an upper bound on the magnitude of the plasma anisotropy parameter since for fixed (p^) (no 
transverse expansion/contraction) ^ cannot be larger than the free-streaming value by causality. 



B. Momentum-space Broadening in a 1+1 Dimensionally Expanding Plasma 

In the previous two subsections we presented details of the limiting cases for 1+1 dimensional plasma evolution: 
1+1 ideal hydrodynamic expansion and 1+1 dimensional free streaming, with the former arising if there is rapid 
thermalization of the plasma and the latter arising if the plasma has no interactions. We would now like to extend 
these models to include the possibility of momentum-space broadening of the plasma partons due to interactions (hard 
and soft). This can be accomplished mathematically by generalizing our expression for ^(r) to 

^(^'-5)= (^^^-1. (3.18) 

In the limit (5 — )■ 0, ^ — > and one recovers the 1+1 hydrodynamical expansion limit and in the limit S 2 one 
recovers the 1+1 dimensional free streaming limit, ^ — > ^f.g. For general S between these limits one obtains the proper 
time dependence of the energy density and temperature by substituting ()3.18|) into the general expression for the 
factorized energy density p.2p to obtain £{t,6) = £oTZ{£,{t, S)). In the limit r 3> tq this gives the following scaling 
relations for the number density, energy density, and hard momentum scale 

71 (t) 
Phard(T) 

£{t) 

Different values of S arise dynamically from the different processes contributing to parton isotropization. Below we 
list the values of S resulting from processes which are relevant during the earliest times after the initial nuclear impact. 



no (-) , (3.19a) 



T 



/To\(l-'5/2)/3 

To (-) , (3.19b) 



T 



r 



1. Collisional Broadening via Elastic 2<-^2 collisions 

In the original version of the bottom up scenario Q, it was shown that, even at early times after the nuclear 
impact, elastic collisions between the liberated partons will cause a broadening of the longitudinal momentum of the 



^ In the original article by Kapusta et. al \3l\ , the calculation of dN^ ' /dydM^d^pT was not presented. 
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particles compared to the non-interacting, free-streaming case. During the first stage of the bottom-up scenario, when 
1 ^ QsT <C as^^, the initial hard gluons have typical momentum of order Qs and occupation number of order l/as- 
Due to the fact that the system is expanding at the speed of light in the longitudinal direction Ng ~ Ql/{asQsT). If 
there were no interactions this expansion would be equivalent to 1-1-1 free streaming and the longitudinal momentum 
PL would scale like l/r. However, when elastic 2 <-> 2 collisions of hard gluons are taken into account [6|, the ratio 
between the longitudinal momentum pL and the typical transverse momentum of a hard particle pT decreases as: 

III oc (Q.r)-2/3 . (3.20) 

Assuming, as before, isotropy at the formation time, tq — Qj^, this implies that for a coUisionally- broadened plasma 
S — 2/3. Note that, as obtained in Ref the derivation of this result makes an implicit assumption that the elastic 
cross-section is screened at long distances by an isotropic real-valued Debye mass. This is not guaranteed in an 
anisotropic plasma as the Debye mass can be become complex due to the chromo-Weibel instability [l3|- However, 
at times short compared to the time scale where plasma instabilities become important we expect the isotropic result 
to hold to good approximation. 



2. Effect of Plasma Instabilities 



Plasma instabilities affect the first stage of bottom-up scenario [I4I. These instabilities are characterized by the 
growing of chromo-electric and -magnetic fields E°- and B"^ . These fields bend the particles and how much bending 
occurs will depend on the amplitude and domain size of the induced chromofields. Currently, the precise parametric 
relations between the amount of plasma anisotropy and amplitude and domain size of the chromofields are not known 
from first principles. There are three possibilities for how the chromo-Weibel instability will affect isotropization of a 
QGP proposed in the literature [4^. l45l[46j: 



(QsT) 



(3.21) 



where 



Ref. [44 1 

45| 



iy = < I Ref. 



(3.22) 



These results correspond to S ~ 1/2, S 



2 Nielsen- Olesen limit, Ref. 46 1 
1/4, and 5 — 1/6, respectively. 



3. Summary and Discussion 



Summarizing, the coefficient S takes on the following values 



2 Free streaming expansion , 
2/3 Collisional-Broadening, Ref. 



1/2 Ref. [44 1 

45| 



1/4 Ref. 
1/6 Nielsen-Olesen limit, Ref. 
Hydrodynamic expansion 



(3.23) 



46 



The exponents in Eq. (|3.23p are direct consequence of the relation between the anisotropy parameter ^ and the 
longitudinal and transverse momentum given in Eq. (|3.ip . The exponent S indicates which kind of broadening we are 
considering. Notice that 6=2 (0) reproduces the behaviour of free streaming (hydrodynamic) expansion. 

In Fig. [3] we sketch the time-dependence of the plasma anisotropy parameter indicating the time scales at which 
the various processes become important. At times shorter than the mean time between successive elastic scatterings, 
'''MFP, the system will undergo 1-1-1 dimensional free streaming with 6 = 2. For times long compared to tmfp but 
short compared to Tinstabiiity the plasma anisotropy will grow with the coUisionally-broadened exponent of 6 = 2/3. 
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FIG. 3: Sketch of the time dependence the plasma anisotropy indicating the various time-scales and processes taking place. 
Here tmfp is the time between elastic collisions (mean-free time) and rinstabiiity is the time at which plasma-instability induced 
soft modes have grown large enough to affect hard particle dynamics. 



Here rinstabmty is the time at which instabiUty-induced soft gauge fields begin to influence the hard-particles' motion. 
When rinstabiiity < T < Tiso the plasma anisotropy grows with the slower exponent of i5 — 1/6... 1/2 due to the 
bending of particle trajectories in the induced soft-field background. At times large compared to rinstabiiity inelastic 
processes are expected to drive the system back to isotropy V\. We note here that for small ^ and realistic couplings 
it has been shown [19] that one cannot ignore the effect of coUisional-broadening of the distribution functions and 
that this may completely eliminate unstable modes from the spectrum. 

Based on such a sketch one could try to construct a detailed model which includes all of the various time scales and 
study the dependence of the process under consideration on each. However, due to the current theoretical uncertainties 
in each of these time scales and their dependences on experimental conditions we choose to use a simpler approach in 
which we will construct two phenomenological models which smoothly interpolate the coefficient S: 

Free streaming interpolating model : 2 > 6 > , 

2 

CoUisionally-broadened interpolating model : - > S > . 

3 

In both models we introduce a transition width, 7~^, which governs the smoothness of the transition from the initial 
value of (5 G {2,2/3} to 5 = at r ~ riso. The free streaming interpolating model will serve as an upper-bound on 
the possible effect of early time momentum-space anisotropics while the collisionally-broadened interpolating model 
should provide a more realistic estimate of the effect due to the lower anisotropics generated. This will help us 
gauge our theoretical uncertainties. Note that by using such a smooth interpolation one can achieve a reasonable 
phenomenological description of the transition from non-equilibrium to equilibrium dynamics which should hopefully 
capture the essence of the physics. In the next section we will give mathematical definitions for these two models. 



C. Space-Time Interpolating Models with Fixed Initial Conditions 



In order to construct our interpolating models, the parameter S should be a function of proper time. To accomplish 
this, we introduce a smeared step function 



A(r,riso,7) 



1 



tanh 



1 



(3.24) 
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FIG. 4: Temporal evolution using our fixed initial condition free-streaming interpolating model {5 = 2) for the energy density 
(left column), hard momentum scale (middle column), and anisotropy parameter (right column) for four different isotropization 
times Tiso £ {1,4,6, 18} tq. The transition width is taken to be (top row) 7 = 2 and (bottom row) 7 — 0.5. To convert to 
physical scales use to ~ 0.3 fm/c for RHIC and to ~ 0.1 fm/c for LHC. 



where 7"^ sets the width of the transition between non-equilibrium and hydrodynamical evolution in units of rjso.* 
In the limit when r <C Tiso, we have A — > and when r ^ riso we have A — > 1. 

Physically, the energy density £ should be continuous as we change from the initial non-equilibrium value of S to 
the final isotropic 6 = value appropriate for ideal hydrodynamic expansion. Once the energy density is specified 
this immediately gives us the time dependence of the hard momentum scale. We find that for general 6 this can be 
accomplished with the following model 



with 7^(0 defined in Eq. 



and for fixed initial conditions 

3A(r)/4 



1 , 



U{t) ^ 
U{t) ^ 



7^((Tiso/To)^-l 

U{t)/U{t^) . 



T 



l-5(l-A(r))/2 



(3.25a) 
(3.25b) 
(3.25c) 



(3.26a) 
(3.26b) 



The power of TZ in hi keeps the energy density continuous at t = Tiso for all 7. In the following subsections we will 
briefiy discuss the two interpolating models we consider in this work. 



Note that compared to Ref. Isoll we have modified our definition of 7 so that it now measures the width in units of Tiso instead of tq. 
This results in time-dependence of modeled quantities not experiencing unphysical "dips" which can occur for large values of 7 in our 
previous interpolating model [48)| . 
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FIG. 5: Temporal evolution using our fixed initial condition coUisionally-broadened interpolating model {S — 2/3) for the 
energy density (left column), hard momentum scale (middle column), and anisotropy parameter (right column) for four different 
isotropization times Tiso £ {1,4,6, 18} tq. The transition width is taken to be (top row) 7 = 2 and (bottom row) 7 = 0.5. To 
convert to physical scales use to ~ 0.3 fm/c for RHIC and tq ~ 0.1 fm/c for LHC. 



1. Free streaming interpolating model 

Using Eq. (|3.25p we can obtain a model which interpolates between early-time 1+1 dimensional longitudinal free 
streaming and late-time 1-1-1 dimensional ideal hydrodynamic expansion by choosing 6 — 2. With this choice and in 
the hmit r <C Tiso, we have A ^ and the system undergoes 1-1-1 dimensional free streaming. When t 3> Tjso then 
A — > 1 and the system is expanding hydrodynamically. In the hmit 7— >oo, A— >0(t — Tjso), the system makes a 
theta function transition from free streaming to hydrodynamical evolution with the energy density being continuous 
during this transition by construction. In Fig. [4] we plot the time-dependence of £, Phaid, and ^ assuming (top) 7 = 2 
and (bottom) 7 = 0.5 for different values of Tiso. As can be seen from this figure for fixed initial conditions during 
the period of free-streaming evolution the system always has a higher effective temperature (phard) than would be 
obtained by a system which undergoes only hydrodynamic expansion from the formation time. As we will show in the 
results section, for fixed initial conditions, this results in a sizable enhancement in high-energy dilepton production. 



2. CoUisionally-broadened interpolating model 

Similarly using Eq. (j3.25p we can obtain a model which interpolates between early-time 1+1 dimensional 
coUisionally-broadened expansion and late-time 1+1 dimensional ideal hydrodynamic expansion by choosing d — 2/3. 
In Fig. Owe plot the time-dependence oi £, piiard, and ^ assuming (top) 7 = 2 and (bottom) 7 ~ 0.5 for different values 
of Tiso. As in the free-streaming interpolating model for fixed initial conditions at early times a coUisionally-broadened 
system always has a higher effective temperature (piiaid) than would be obtained by a system which undergoes only 
hydrodynamic expansion from the formation time. As we will show in the results section, for fixed initial conditions, 
this results in an enhancement in high-energy dilepton production; however, compared to the free-streaming case the 
effect is reduced due to the lower effective temperatures obtained by the coUisionally-broadened plasma. We also note 
that in the case of coUisionally-broadened expansion the magnitude of ^ is significantly reduced as compared to the 
free-streaming case. As can be seen from the rightmost panel of Fig. [5] even if one assumes a large isotropization time, 
Tiso = 18 To, the amount of momentum space anisotropy generated is small with ^max ~ 2.5 for 7 = 2 and ^max ~ 1-5 
for 7 = 0.5. 
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FIG. 6: Percentage entropy generation using fixed initial condition interpolating models (|3.25p with 5 = 2 and 5 = 2/3. 
Horizontal lines show 10% and 20% entropy generation bounds. 



Interpolating Model 


RHIC ~ 10% 


RHIC - 20% 


LHC - 10% 


LHC - 20% 


Free-Streaming (5 — 2) 


Tiso < 0.8 fm/c 


Tiso < 1.2 fm/c 


Tiso < 0.26 fm/c 


Tiso < 0.4 fm/c 


CoUisionally-Broadened (5 = 2/3) 


Tiso < 5 fm/c 


Tiso < 18 fm/c 


Tiso < 1.6 fm/c 


Tiso < 6.2 fm/c 



TABLE I: Bounds on riso imposed by requiring either a 10% or 20% bound on percentage entropy (particle number) generation 
from our fixed initial condition interpolating models. To convert to physical scales we have used tq = 0.3 fm/c for RHIC and 
TO = 0.1 fm/c for LHC. 



D. Space-Time Interpolating Models with Fixed Final Multiplicity 



In the previous subsection we constructed models which allow one to interpolate between an initially non-equilibrium 
plasma to an isotropic equilibrium one assuming that the initial conditions are held fixed. One problem with this 
procedure is that given fixed initial conditions these interpolating models will result in generation of particle number 
during the transition from S G {2,2/3} to zero. 

One can derive an expression for the amount by which the number density is increased by starting from the general 
expression for the particle number density n(T)/no = (phard/2o) (1 + £,{t))^^^^ and using the expression for phard 
derived in the previous section p.25cp to obtain 



n{T) U{t) 



Taking the limit t ^ Tjso we obtain 



lim 



Mr) 



ISO / ' ISO 



5/2-1 



7^ (Wro) -1 



3/4 



-r»Tiao no T \ Tq 

Translating this into a statement about the entropy generation using S{t) — Tn{T) gives 



r>Tiao bo 



To 



5/2 



n 



{{Tiso/ToY 



1 



3/4 



(3.27) 



(3.28) 



(3.29) 



When either S ^ or riso tq, AS" = (5'finai — So) /So goes to zero and there is no entropy generation; however, 
entropy generation increases monotonically with 5. In the limit of large Tiso/ To we find 



lim ^ 

iBo-»oo Do 



To 



- 1 



(3.30) 
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FIG. 7: Temporal evolution using our fixed final multiplicity interpolating models for the energy density (left column), hard 
momentum scale (middle column), and anisotropy parameter (right column) for four different isotropization times riso £ 
{1, 4, 6, 18} To. Top row is the free-streaming interpolation model (5 = 2) and bottom row is the coUisionally-broadening 
interpolation model {5 — 2/3). To convert to physical scales use to ~ 0.3 fm/c for RHIC and tq ~ 0.1 fm/c for LHC. 



Again we see that in the limit that either (5 — > or Tiso tq then there is no entropy generation. 

The requirement of bomided entropy generation can be used to constrain non-equilibrium models of the QGP [4^ . 
In Fig. [HI we plot the entropy generation (particle number generation) resulting from our models using Eq. (j3.29p for 
6 £ {2, 2/3} along with bounds at 10% and 20%. In the free-streaming interpolating model {6 — 2) with fixed initial 
conditions requiring that the percentage entropy generation be less than each of these bounds requires Tjso < 2.6 tq 
for the 10% bound and Tiso < 4ro for the 20% bound. In the coUisionally-broadened interpolating model {6 = 2/3) 
with fixed initial conditions we obtain similarly Tiso < 17 tq for the 10% bound and Tiso < 62 tq for the 20% bound. 
We summarize our results in Table HI As can be seen from Table |T] requiring the listed bounds on entropy generation 
the values of Tiso allowed in our free-streaming interpolating model become highly constrained. However, in the case 
of the coUisionally-broadened interpolating model the upper-bounds imposed on Tiso are much larger due to the much 
lower entropy generation required to transition from coUisionally-broadened evolution to hydrodynamic evolution. 

One problem with our fixed initial condition family of models is that due to the fact that they generate additional 
particles the multiplicity of final particles is not independent of the assumed value of Tjso- Because most of the 
experimental results for dilepton spectra are binned with respect to a fixed final multiplicity this means that we 
should also construct models which always result in a fixed final number density. In the following subsection we will 
show how this can be accomplished. 



1. Enforcing Fixed Final Multiphcity 



We will now construct interpolating models which have a fixed final entropy (multiplicity). In order to accomplish 
this the initial conditions will have to vary as a function of the assumed isotropization time. We will show that, as a 
result, for finite Tiso one must lower the initial "temperature" in both the free-streaming and coUisionally-broadened 
interpolating models. To accomplish this requires only a small modification to the definition of U in Eq. (j3.26p : 



lim U{t) 

+ 



n 



3/4 



To 



(3.31a) 
(3.31b) 
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As a consequence of this modification, the initial energy density and hence initial "temperature" will depend on the 
assumed value for Tiso- There is no modification required for ^. We demonstrate this in Fig. [7] where we plot the 
time-dependence of Phmd, and ^ for 7 = 2 and (top) S — 2 and (bottom) 6 — 2/3. 

In the remainder of this work we present our final results for dilcpton yields using both approaches, i.e., fixed 
initial conditions using Eqs. ((3?25)) with ((3?26)) or fixed final multiphcity through Eq. ((3?25l) with ((3?3T|) . We mention 
that in both cases, dilepton production is affected in the presence of anisotropics in momentum-space, however, one 
anticipates that the effect will be larger when the initial conditions are held fixed due to the larger particle number 
generation. We will come back to this issue in the conclusions and discussion. 

IV. RESULTS 

In this section we will present expected e+e~ yields resulting from a central Au-Au collision at RHIC full beam 
energy, -^=200 GeV and from a Pb-Pb collision at LHC full beam energy, -^5=5. 5 TeV. In all figures in this section 
we will present the prediction for RHIC energies in the left panel and LHC energies in the right panel. 

Before presenting our results we first explain the setup, numerical techniques used, and parameters chosen for our 
calculations. Because the differential dilepton rate di?'*' /d'^P given in Eq. (|2.8p is independent of the assumed 
space-time model. We first evaluate it numerically using double-exponential integration with a target precision of 
10~^. The result for the rate was then tabulated on a uniformly-spaced 4-dimensional grid in M, pT, y, and logj^g C 
: M/phard,PT/Phard & {0.1,25}, y G {—3,3}, log^Q ^ G {—6,4}. This table was then used to build a four-dimensional 
interpolating function which was valid at continuous values of these four variables. We then boost this rate from 
the local reference frame to center-of-mass frame and evaluate the remaining integrations over space-time (r and 
77) and transverse momentum or invariant mass appearing in Eqs. (j2.1ip using quasi-Monte Carlo integration with 
r € {ro,T^}, r/ G {—2.5,2.5} and, depending on the case, restrict the integration to any cuts specified in M or px- 
Our final integration time, Tf, is set by solving numerically for the point in time at which the temperature in our 
interpolating model is equal to the critical temperature, i.e. Phaid('''/) = 7c. We will assume that when the system 
reaches Tc, all medium emission stops. We are not taking into account the emission from the mixed/hadronic phase 
at late times since the kinematic regime we study (high M and pt) is dominated by early-time high-energy dilepton 
emission [sol. l33l|. 

For RHIC energies we take an initial temperature To— 370 MeV, at a formation time of ro= 0.26 fm/c, and use 
Rt— 6.98 fm. For LHC energies, we use to= 0.088 fm/c, Tq= 845 MeV and Rt— 7.1 fm. In both cases, the critical 
temperature Tc is taken as 160 MeV and the spectra are calculated at midrapidity region y = 0. Any cuts in 
transverse momentum or invariant mass will be indicated along with results. Note that the precise numerical value 
of the parameters above were chosen solely in order to facilitate straightforward comparisons with previous works 
[SS*! from which we have obtained predictions for Drell Yan, heavy quark, jet-fragmentation, and jct-thcrmal dilcpton 
yields. 

Finally we note that below we will use i^-factors to adjust for next-to-leading order corrections to the dilepton 
rate. These /C-factors are determined by computing the ratio of the next-to-leading order prediction of [H, [sl] with 
our leading order prediction in the case of ideal-hydrodynamic expansion. We therefore assume that the if-factors 
are independent of the assumed thermalization time. This is an approximation which, in the future, one would like 
to relax by computing the full next-to-leading order dilepton rate in the presence of momentum-space anisotropics. 

A. Dilepton production with fixed initial conditions 

We now present the results of the dilepton production assuming the time dependence of the energy density, the 
hard momentum scale and the anisotropy parameter are given by Eqns. (|3.25p and (|3.26p with 6 G {2, 2/3}. 

1. Free streaming interpolating model 

In Fig. [HI we show our predicted dilepton mass spectrum for RHIC and LHC energies assuming the time dependence 
of the energy density, the hard momentum scale and the anisotropy parameter are given by Eqns. p.25p and (j3.26p with 
d — 2. This corresponds to our free-streaming interpolating model. For us this model will serve as an upper-bound on 
the possible effect of momentum-space anisotropics on dilepton yields. From Fig. [5]we see that for both RHIC or LHC 
energies, there is a significant enhancement of up to one order of magnitude in the medium dilepton yield when we vary 
the isotropization time from tq to 2 fm/c. This enhancement is due to the fact that in 1-1-1 dimensional free streaming, 
the system preserves more transverse momentum as can be seen from Fig. [3) For fixed initial conditions this results 
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FIG. 8: Free-streaming interpolating model dilepton yields as a function of invariant mass in central Au-j-Au collisions at RHIC 
(left) and Pb-I-Pb at the LHC (right), with a cut pt > 4 (8) GeV and rapidity y=0. For medium dileptons we use 7=2 and riso 
is taken to be either 0.26 (0.088) fm/c or 2 fm/c for RHIC (LHC) energies with fixed initial conditions. A A'-factor of 1.5 was 
applied to account for NLO corrections. Dilepton yields from Drell Yan, Heavy Quarks, Jet-Thermal and Jet-Fragmentation 
were obtained from Ref. [331. 
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FIG. 9: Free-streaming interpolating model dilepton yields as a function of transverse momentum in central Au+Au collisions 
at RHIC (left) and Pb+Pb at the LHC (right), with a cut 0.5 < M < 1 GeV and rapidity y—Q. For medium dileptons we 
use 7=2 and riso is taken to be either 0.26 (0.088) fm/c or 2 fm/c for RHIC (LHC) energies with fixed initial conditions. A 
/S'-factor of 6 was applied to account for NLO corrections. Dilepton yields from Drell Yan, Jet- Thermal and Jet-Fragmentation 
were obtained from Ref. [38j . 



in a larger effective temperature than would be obtained if the system underwent locally-isotropic (hydrodynamical) 
expansion throughout its evolution. 

Nevertheless, as Fig. [5]shows, as a function of invariant mass, the other contributions to high-energy dilepton yields 
(Drell- Yan, jet-thermal, and jet-fragmentation) are all of the same order of magnitude as the medium contribution. 
This coupled with the large background coming from semileptonic heavy quarks decays would make it extremely 
difhcult for experimentalists to extract a clean medium dilepton signal from the invariant mass spectrum. For this 
reason it does not look very promising to determine plasma initial conditions from the dilepton invariant mass 
spectrum. For this reason we will not present our predictions for the invariant mass spectrum for the intermediate 
models detailed below and only return to the invariant mass spectrum at the end of this section for completeness. 

The good news is, however, that as a function of transverse momentum, see Fig. ^ the production of medium 
dileptons is expected to dominate other production mechanisms for pr < i (6) GeV in the case of RHIC (LHC). 
In addition to this we see that for the free-streaming interpolating model that there is a significant enhancement of 
medium dileptons for both RHIC and LHC energies. 



18 



4Eh Drell Yan 




123456 2 4 6 8 10 

Pt [GeV] Pt [GeV] 

FIG. 10; CoUisionally-broadened interpolating model dilepton yields including coUisional broadening as a function of transverse 
momentum in central Au+Au collisions at RHIC (left) and Pb+Pb at the LHC (right), with a cut 0.5 < M < 1 GeV and 
rapidity y—Q. For medium dileptons we use 7=2 and riso is taken to be either 0.26 (0.088) fm/c or 2 fm/c for RHIC (LHC) 
energies with fixed initial conditions. A A"-factor of 6 was applied to account for NLO corrections. Dilepton yields from Drell 
Yan, Jet-Thermal and Jet-Fragmentation were obtained from Ref. [S^. 
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FIG. 11: Dilepton enhancement, 0, as defined in Eq. (|4.1|l resulting from our coUisionally-broadened interpolating model 
((5 = 2/3) with fixed initial conditions and tiso = 2 fm/c. Result for RHIC energies is shown on left and for LHC energies 
on right. The invariant mass cut used was 0.5 < M < 1 GeV and rapidity y=0. Lines show expected pre-equilibrium 
dilepton enhancements for different values of the transition width 7 corresponding to sharp or smooth transitions between 
pre-equilibrium and equilibrium behavior. 



In order to quantify the effect of time-dependent pre-equilibrium emissions we define the "dilepton enhancement" , 
0('''iso), as the ratio of the dilepton yield obtained with an isotropization time of Tjso to that obtained from an 
instantaneously thermalized plasma undergoing only 1-1-1 hydrodynamical expansion, ie. Tiso — tq. 

^(^. ) ^ ( dN^^^JnsoA I f diV^^"'(riso = ro) 
y dydp'^ j I \ dydp^ 

Using this criterion wc find for the free streaming interpolating model with fixed initial conditions the dilepton 
enhancement at Tiso = 2 fm/c can be as large as 10. However, as mentioned above we expect that the actual 
enhancement will be lower due to the fact that parton interactions such as coUisional-broadening will modify the 
free-streaming ^ = t^/tq — 1 to something growing slower in proper time bringing the system closer to equilibrized 
expansion. In addition, as we will discuss below when using fixed initial conditions and 5 = 2 there is significant 
entropy generation which, when properly normalized to fixed final multiplicity, results in reduced (j). Therefore, we 
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FIG. 12: Free streaming interpolating model dilepton yields as a function of transverse momentum in central Au+Au collisions 
at RHIC (left) and Pb+Pb at the LHC (right), with a cut 0.5 < M < 1 GeV and rapidity y—Q. For medium dileptons we 
use 7=2 and riso is taken to be either 0.26 (0.088) fm/c or 2 fm/c for RHIC (LHC) energies and fixed final multiplicity. A 
/('-factor of 6 was applied to account for NLO corrections and rapidity y=0. Dilepton yields from Drell Yan, Jet- Thermal and 
Jet-Fragmentation were obtained from Ref. [3^ . 

expect ^ 10 obtained from the free streaming interpolation model with fixed initial conditions to be an upper-bound 
on the effect of pre-equilibrium emissions. Some of our results fixing initial conditions are related with recent work 
on dilepton production from a viscous QGP [1^. 

2. Collisionally-broadened interpolating model 

In Fig. 1101 we show our predicted dilepton transverse momentum spectrum for RHIC and LHC energies assuming 
the time dependence of the energy density, the hard momentum scale and the anisotropy parameter are given by 
Eqns. (|3.25p and (|3.26p with (5 = 2/3. This corresponds to our collisionally-broadened interpolating model with fixed 
initial conditions. From Fig. (TU] we see that for both RHIC or LHC energies, there is dilepton enhancement in the 
kinematic range shown; however, compared to the free streaming case the enhancement is reduced. This is due to 
the fact that the collisionally-broadened interpolating model is always closer to locally-isotropic expansion than the 
free-streaming (^ = 2) model, see Figs. |4] and [3 

In Fig. [TT] we show the dilepton enhancement, 0, as function of transverse momentum for Tjso = 2 fm/c at (left) 
RHIC energies (right) LHC energies. The invariant mass cut is the same as in Fig. [10] (0.5 < Af < 1 GeV). As 
can be seen from Fig. [TT] using fixed initial conditions there is a rapid increase in (f) between 1 and 3 GeV at RHIC 
energies and 1 and 4 GeV at LHC energies. The precise value of the enhancement depends on the assumed width 7"^ 
and in Fig. [TT] we show for 7"^ S {0.1, 0.5, 2}. As can be seen from this figure both sharp and smooth transitions 
from early-time collisionally-broadened expansion to ideal hydrodynamic expansion result in a 40-70% enhancement 
of medium dilepton yields at RHIC energies and 60-100% at LHC energies. We will return to this in the results 
summary at the end of this section. 

B. Dilepton production with fixed final multiplicity 

We now present the results of the dilepton production assuming the time dependence of the energy density, the 
hard momentum scale and the anisotropy parameter are given by Eqns. (|3.25p and (|3.3ip with 5 G {2, 2/3}. 

1. Free streaming interpolating model 

In Fig. I12[ we show our predicted dilepton transverse momentum spectrum for RHIC and LHC energies assuming 
the time dependence of the energy density, the hard momentum scale and the anisotropy parameter are given by 
Eqns. (|3.25p and (|3.3ip with 5 = 2. This corresponds to our free-streaming interpolating model with fixed final 
multiplicity. From Fig. [T^] we see that for both RHIC or LHC energies, there is dilepton enhancement in the 



20 



> 



T3 



O 



-9 
-10 
-11 





Heavy Quarks 


-B- 


Drell Yan 




Jet-Thermal 




Jet Fragmentation 


▲ 


Medium (riso = 0.26 fm/c) 


■ 


Medium (rjso = 2 fm/c) 



3 4 
M [GeV] 




> 
o 
O 



■a 
>. 



o 





Heavy Quarks 




Drell Yan 




Jet-Medium 




-A- Medium (rbs = 0.088 fm/c) 


' 1 

^ 1 , 


-■- Medium (Tjso = 2 fm/c) 






























6 

M [GeV] 



10 



FIG. 13: CoUisionally-broadened interpolating model dilepton yields as a function of invariant mass in central Au+Au collisions 
at RHIC (left) and Pb+Pb at the LHC (right), with a cut pr > 4 (8) GeV and and rapidity j/=0. For medium dileptons we 
use 7=2 and Tiso is taken to be either 0.26 (0.088) fm/c or 2 fm/c for RHIC (LHC) energies and fixed final multiplicity. A 
Js'-factor of 1.5 was applied to account for NLO corrections. Dilepton yields from Drell Yan, Heavy Quarks, Jet-Thermal and 
Jet-Fragmentation were obtained from Ref. [3^ . 
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FIG. 14: Collisionally-broadened interpolating model dilepton yields as a function of transverse momentum in central Au+Au 
collisions at RHIC (left) and Pb-|-Pb at the LHC (right), with a cut 0.5 < M < 1 GeV and rapidity y=0. For medium 
dileptons we use 7=2 and Tiso is taken to be either 0.26 (0.088) fm/c or 2 fm/c for RHIC (LHC) energies and fixed final 
multiplicity. A A'-factor of 6 was applied to account for NLO corrections. Dilepton yields from Drell Yan, Jet-Thermal and 
Jet-Fragmentation were obtained from Ref. [3^ . 



kinematic range shown; however, when fixing on final multiplicity the effect of a free-streaming pre-equilibrium phase 
is reduced. In fact, for small and large pT the free-streaming interpolating model with fixed final multiplicities predicts 
a suppression of dileptons. This is due to the fact that in order to maintain fixed final multiplicity for Tjso = 2 fm/c 
the free-streaming model initial energy density has to be reduced by ~ 50% (see top row of Fig. [7]). 



2. Collisionally-broadened interpolating model 



In Figs. 1131 and [HI we show our predicted dilepton invariant mass and transverse momentum spectrum for RHIC 
and LHC energies assuming the time dependence of the energy density, the hard momentum scale and the anisotropy 
parameter are given by Eqns. (I3.25P and p.3ip with S = 2/3. This corresponds to our collisionally-broadened 
interpolating model with fixed final multiplicity. In Fig. [T5| we show the dilepton enhancement, 0, as function of 
transverse momentum for Tiso = 2 fm/c at (left) RHIC energies (right) LHC energies. The invariant mass cut is 
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FIG. 15: Dilepton enhancement, 0, as defined in Eq. (|4.1|l resulting from our coUisionally-broadened interpolating model 
{5 = 2/3) with fixed final multiplicity and riso = 2 fm/c. Result for RHIC energies is shown on left and for LHC energies 
on right. The invariant mass cut used was 0.5 < M < 1 GeV and and rapidity y=Q. Lines show expected pre-equilibrium 
dilepton enhancements for different values of the transition width 7 corresponding to sharp or smooth transitions between 
pre-equilibrium and equilibrium behavior. 



the same as in Fig. [T3] (0.5 < M < 1 GeV). As can be seen from Fig. [T5] similar to the case of fixed initial 
conditions there is a rapid increase in </> between 1 and 3 GeV at RHIC energies and 1 and 4 GeV at LHC energies. 
However, compared to the case of the coUisionally-broadened interpolating model with fixed initial condition (Fig. fTTj) 
the maximum enhancement is reduced slightly and we see a more pronounced peak in (/> as a function of transverse 
momentum appearing. As can be seen from this figure both sharp and smooth transitions from early-time coUisionally- 
broadened expansion to ideal hydrodynamic expansion result in a 20-40% enhancement of medium dilepton yields at 
RHIC energies, and 30-50% at LHC energies. 



C. Summary of Results 



Based on the figures presented in the previous subsections we see that the best opportunity for measuring information 
about plasma initial conditions is from the M < 2 GeV dilepton transverse momentum spectra between 1 < < 
6 GeV at RHIC and 2 < pT < 8 GeV at LHC. This is due to the fact that medium dilepton yields dominate other 
mechanisms in that kinematic range and hence give the cleanest possible information about plasma initial conditions. 
In all cases shown above dilepton production is enhanced by pre-equilibrium emissions with the largest enhancements 
occurring when assuming fixed initial conditions and the free-streaming interpolating model. As we have mentioned 
above this model sets the upper-bound for the expected dilepton enhancement. Our most physically realistic model is 
the coUisionally-broadened interpolating model with fixed final multiplicity so we will use it for our final predictions 
of expected dilepton enhancement. For this model, as can be seen from Fig. 1151 assuming Tjso — 2 fm/c we find a 
20-40% enhancement in dilepton yields at RHIC and 30-50% at LHC. 

In addition we can calculate the dilepton enhancement for different assumed values for Tiso. This is shown for RHIC 
energies (left) and LHC energies (right) in Fig. [1^] where we have fixed 7 = 2 and varied Tiso to see the effect of 
varying the assumed isotropization time. As can be seen from this figure the effect of reducing Tjso is to shift the peak 
in (j) to larger pT while at the same time reducing the overall amplitude of the peak. This feature seems generic at 
both RHIC and LHC energies. Therefore, in order to see the difference between an instantaneously thermalized QGP 
with Tiso = To and one with a later thermalization time requires determining the medium dilepton spectra between 
1 < Pt < 6 GeV at RHIC and 2 < < 8 GeV at LHC with high precision so that one could measure the less than 
50% variation resulting from pre-equilibrium emissions. 

Finally, we point out that in Fig. [11] we have chosen an invariant mass cut of 0.5 < M < 1 GeV. Since our model 
predicts the full yields versus Af and pr it is possible to take other cuts (invariant mass and/or transverse momentum). 
This could be coupled with fits to experimental data, allowing one to fix Tjso and 7 via a "multiresolution" analysis. 
To demonstrate the dependence of (j) on the mass cut in Fig. [17] we show the dilepton enhancement, (j), using a mass 
cut of 1 < M < 2 GeV. As can be seen from this Figure the qualitative features of our model's predictions are similar 
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FIG. 16: Dilepton enhancement, 0, as defined in Eq. (|4.1|l resulting from our collisionally-broadened interpolating model 
(5 = 'i/'i) with fixed final multiplicity and 7 = 2. Result for RHIC energies is shown on left and for LHC energies on right. 
The invariant mass cut used was 0.5 < M < 1 GeV and and rapidity y=G. Lines show expected pre-equilibrium dilepton 
enhancements for different values of the assumed plasma isotropization time, Tiso. 
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FIG. 17: Dilepton enhancement, 0, as defined in Eq. (|4.ip resulting from our collisionally-broadened interpolating model 
[5 — 2/3) with fixed final multiplicity and 7 = 2. Result for RHIC energies is shown on left and for LHC energies on right. 
The invariant mass cut used was 1 < M < 2 GeV and and rapidity j/=0. Lines show expected pre-equilibrium dilepton 
enhancements for different values of the assumed plasma isotropization time, Tiso- 



to the lower mass cut presented in Fig. 1161 however, for this mass cut we see that there is a stronger suppression of 
dilepton production at low and high invariant masses if there is late thermalization, Tjso > 2 fm/c. Such features can 
be used to constrain the model further when confronted with experimental data. 



V. CONCLUSIONS 



In this work we have presented models which allow one to smoothly interpolate between early-time non-equilibrium 
1+1 dimensional expansion to late-time isotropic equilibrium 1+1 dimensional hydrodynamic expansion. To accom- 
plish this we introduced simple interpolating models with two parameters: Tjso, which is the time at which the system 
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begins to expand hydrodynamically and 7 which sets the width of the transition. Using these models we integrated 
the leading order rate for dilepton production in an anisotropic plasma over our modeled space-time evolution. Based 
on our numerical results for the variation of dilepton yields with the assumed values of riso we find that the best 
opportunity to determine information about the plasma isotropization time is by analyzing the high transverse mo- 
mentum {1 < pt < 6 GeV at RHIC and 2 < pr < 8 GeV at LHC) dilepton spectra using relatively low pair invariant 
mass cuts {M < 2 GeV). Based on these px spectra we introduced the "dilepton enhancement" factor (/)(Tiso) which 
measures the ratio of yields obtained from a plasma which isotropizes at riso to one which isotropizes at the formation 
time, tq. 

We showed that for our most extreme model, the free-streaming interpolating model {6 = 2) with fixed initial 
conditions, that the resulting enhancement can be as large as 10; however, this extreme model probably overestimates 
the amount of anisotropy in the plasma. Additionally, this model results in a large amount of entropy generation 
during the transition from the free-streaming asymptotic behavior to hydro r~^/^ asymptotic behavior. As 
we discussed this greatly constrains the maximum isotropization times Tjso which are consistent with experimental 
indications of low (10-20%) entropy generation. 

In order to construct a more realistic model we then included collisional-broadening of the initial pre-equilibrium 
parton distribution functions {S = 2/3). In this more realistic model there is much less entropy generation and the 
system is always closer to ideal 1-1-1 hydrodynamic expansion than in the free-streaming interpolating model. As 
a result the dilepton enhancement due to prc-equilibrium emissions is lower than the free-streaming case. We find 
that when fixing final multiplicity at RHIC energies there is a 20-40% enhancement in the high-transverse momentum 
dileptons and at LHC energies it is 30-50% when one assumes an isotropization time of Tiso — 2 fm/c. The amplitude of 
the enhancement and position of the peak in the enhancement function, (j), varies with the assumed value of Tiso which, 
given sufficiently precise data, would provide a way to determine the plasma isotropization time experimentally. We 
presented our predictions for the dilepton enhancement, 0, as a function of riso for two different invariant mass cuts, 
demonstrating that our model can be constrained by a multiresolution analysis which should give higher statistics 
and further constrain the two model parameters at our disposal. 

One shortcoming of this work is that we haven't included NLO corrections to dilepton production from an anisotropic 
QGP. At low invariant mass these corrections would become important. As a next step one must undertake a 
calculation of the rate for dilepton pair production at NLO in an anisotropic plasma. This is complicated by the 
presence of plasma instabilities which render some expressions like (AA) correlators formally divergent and hence 
analytically meaningless. However, when combined with numerical solution of the long-time behavior of a plasma 
subject to the chromo-Weibel instability it may be possible to extract finite correlators [i^. This is a daunting but 
doable task. Absent such a calculation, phenomenologically speaking it is probably a very good approximation to 
simply take existing NLO calculations and apply the enhancement function (p as calculated at LO. We leave this for 
future work. 

Another uncertainty comes from our implicit assumption of chemical equilibrium. If the system is not in chemical 
equilibrium (too many gluons and/or too few quarks) early time quark chemical potentials, or fugacities, will affect 
the production of lepton pairs (3^. [33^. However, to leading order the quark and gluon fugacities will cancel between 
numerator and denominator in the dilepton enhancement, (j) (33j . We, therefore, expect that to good approximation 
one can factorize the effects of momentum space anisotropies and chemical non-equilibrium. 

We note in closing that the interpolating model presented here has application beyond the realm of computing 
dilepton yields. In fact, such a model can be used to assess the phenomenological consequences of momentum-space 
anisotropies in other possible observables which are sensitive to early-time stages of the QGP, e.g., photon production 
(50| . heavy-quark transport, jet-medium induced electromagnetic radiation, etc. 
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